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Evidence of finite sQGP formation time at RHIC 
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Abstract. We demonstrate that the existence of a finite formation time 
of strongly interacting plasma in nuclear collisions at RHIC is an inavitable 
conclusion from recent experimental data. The most striking feature of the 
experimental data - an absense of absorption of high transverse momentum 
pions in the rection plane direction for mid-pcriphcral collisions - points to 
the presence of a surface zone with no absorption and strong suppression in 
the inner core. A natural interpretation of such a zone should be the plasma 
formation time T~ 2-3 fm/c. Wc discuss constraints induced by finite forma- 
tion time on some physical observables at RHIC. Nuclear modification factor, 
azimuthal assymetry, di-jet correlations can quantitatively be described by par- 
ticle production in the early stage of the nuclear collision. A possible impact 
on azimuthal anisotropy at lower hardon momenta, on interpretation of the 
non-photonic electrons and J/t/> data is also considered. 

Keywords: jet absorption, corona effect, Woods-Saxon, quark-gluon plasma, 
reaction plane, formation time 
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1. Introduction 

Nowadays, the physics of relativistic nuclear collisions at RHIC goes from the phase 
of discoveries [Q to the stage of understanding properties of a new state of nuclear 
matter: the quark-gluon plasma, QGP. Naive expectations that QGP to be a gas 
of weakly interacting quark and gluons was washed out by very first experimantal 
data. Collective phenomena, such as radial and elliptic flow, reveal "hydro" prop- 
erties of the QGP, which behave "like a good liquid rather than a dilute gas of 
quasiparticles" [E]. The idea was introduced, that plasma at RHIC should be in 
strongly coupled regime, sQGP. Another phenomenon seen at RHIC is jet quench- 
ing. Even very hard jets are expected to lose some energy in the formed medium. 
The quenching is decribed by a nuclear modification factor, Raa, which is defind as 
the observed number of jets normalized to the expected number of jets from super- 
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position of individual nucleon-nucleon collisions. Very first RH1C data have indeed 
shown significant particle suppression at high px [ 1 ■ Recently, the PHENIX col- 
laboration showed a preliminary result for the dependence of Raa on the azimuthal 
angle 4> relative to the reaction plane in Au+Au collisions at v / saT/v=200 GeV 
at RHIC [El EJ Eli Fig. ^ The most interesting feature of these data is that, at 
event centrality class 50-60%, for transverse momenta above 4 GeV/c, in-plane Raa 
equals to one within the errors. This implies no absorption at all for high pt pions. 
For such high momenta, Cronin effect [E] is negligible and can not bring Raa close 
to one. At the same time, a significant particle absorption is seen in out-of-plane. 
Apparently, parton energy loss calculations can't describe this "punch through" 
feature of in-plane Raa=^, e.g. see Fig. 5 in Ref. [0]. In our resent papers [[7\ El 
we offer an alternative explanation of this and certain other features of the data. 
We introduce a new term plasma formation time, which could be considered as a 
transision time for produced matter to reach liquid or strongly interacting phase. 
During this time, soft and hard jets can escape from interaction zone with very 
little or no energy loss. Only jets originating on the surface can reach the detector. 
It forms kind of corona around nuclear overlap region with thickness significantly 
larger than nuclear diffuscness zone. In this paper we describe consequences of 
sQGP formation time on some physical observables. 

2. Where formation time comes from? 

Let me briefly describe the arguments were used in Ref. [ [7j. We construct a 
simple model using the Monte Carlo simulation of nucleus-nucleus collisions based 
on the Glauber approach and discussed in detail in [El- A Woods-Saxon density 
distribution was used in all cases. We restrict ourselves to describe the data of high 
Pt pions with transverse momentum above 4 GeV/c, where Raa does not depend 
on pt- We assume that all high px pions are produced by parton fragmentation and 
that the number of partons is proportional to Nbinary- If there is no absorption, 
Raa=1 in all directions and the shape of the event is isotropic. To explain the 
experimentally observed feature of in-plane Raa—^, we investigate the role of purely 
geometrical factors. In our model, jets which have to travel through the medium 
at some direction from thier production point to the surface less than a distance L 
will leave the interaction zone unmodified. Jets originating in the core region deeper 
than L suffer significant energy loss and are completely absorbed. We can expalin 
the picture in other words: let for each jet to be free streaming at any direction 
with speed of light, c, for some time T. After that time we stop all jets which 
still left within the envelop of Woods-Saxon radia of overlapping nuclei. The whole 
picture looks like a pure corona jet production, but we allow this corona region to 
be larger than a Woods-Saxon type skin. We can calculate Raa as a ratio of the 
seen collisions, N co u, to the total Nbinary at particular centrality class. The cut-off 
parameter L = T * c should be on the order of the size of the in-plane interaction 
zone at 50-60% centality. Parameter L was found to be L=2.3±0.6 fm to get 
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Fig. 1. PHENIX preliminary Raa results for ir° at momenta 5-6 GeV/c versus 
angle <j> relative to the reaction plane for different centralities in Au+Au collisions 
at RHIC [laj- Points are experimental data, thin lines show systematic errors from 
the reaction plain resolution, vertical bars in the middle show avaraged over the 
reaction plane Raa value and it's error. Black horizontal bars are our predictions. 



i?™ A =0.9±0.1 for 50-55% centrality. This is close to what is seen experimentally 
and leaves some room for Cronin enhancement [ |S] , if any. The results of our 
calculations together with experimental data for Au+Au collisions at 200 GeV are 
shown in Fig. ^ and Fig- El In Fig- 13 we also plot our results for Raa at high pt 
in Cu+Cu collisions at 200 GeV using the same L=2.3 fm. There is an axcellent 
agreement with all experimental data. 

We investigated the sensitivity of our result to various assumptions, like: if 
the thickness of the material integrated over the path length is a critical cut-off 
parameter, the centrality dependence becomes very strong in this case and can not 
describe the data. Another premise tested uses the number of participant nucleons, 
Npart, instead of N co u. In this case the centrality dependence is weaker than in the 
experimental data. To take into account obvious fluctuations in the jet absorption 
we smooth the cutoff edge L by introdusing its diffuseness D of 0.1-2 fm. The results 
do not change significantly in this parameter range. A maximum deviation of 5% 
was obtained for D=2 fm, only in very peripheral collisions. The default value was 
choosen to be D=l fm. 
Within our model we can also describe Raa—0-^G for neutral high px pions in 
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Au+Au collisions at v / saT/v=62.4 GeV [0] in most central events as measured by 
PHENIX, but with even larger value T=3.5 fm/c. 




Fig. 2. Raa for Au+Au (dashed curve) and Cu+Cu (solid curve) collisons versus 
the number of participant nucleons, N part . The cicles are experimental 7r° data 
for Au+Au collisions integrated for j>t>4 GeV/c [113). The triangles are data for 
Cu+Cu collisions of pr>7 GeV/c [E|- Only statistical errors are shown. 



3. Discussion and Conclusions 



What could be the physical interpretation of the cutoff T = L/cl Our guess is 
that it is natural to assign this parameter T to a plasma formation time, or, at 
least, the time when parton energy loss actually starts. This gives a simple and 
elegant interpretation of the effect: particles produced close to the surface of the 
collision zone have time, about 2-3 fm/c, to escape. After that time a dense and 
strongly interacting matter, sQGP, is formed and high px partons suffer from a 
significant energy loss. This formation time puts significant constraints on many 
"explained" and not explained yet observables at RHIC. Let us go through the list 
of such observables and some physical consequences. 

1) Formation time serves as a corona for particle production at the early stage 
of the collision. So, we have to split particle production into two sources, from 
the core and from the corona. Most of soft particle are coming from core, all hard 
particles are produced in corona. There should be the later stage hadron gas corona 
too, but it might have very little influence on high px particles. 

2) The core region is very opaque, all jets are stopped in the core. Indeed, there 
are some theoretical estimations that jets would not make it beyond 1/3 or 1/2 fin [ 

ma Eg. 

3) All pions with momentum above 3-4 GeV/c, all baryons - above 5-6 GeV/c are 
coming purely from corona, this is why Raa is constant at high pt [ EH • It also 
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Fig. 3. Calculated cllipticity parameter v^ for Au+Au collisons, solid line, versus 
the number of participant nucleons, N part . Data for ir° with error bars are: cicles 
for 4.59 GeV/c, squares for 5-7 GeV/c. PHENIX preliminary data [HU 0. 



means that the information about sQGP properties lays solely on low momenta 
particles. 

4) The experimental variable, v%, describes the elliptic shape of the event in the 
nucleus-nucleus collision in azimuthal angle. We estimate V2 for high px parti- 
cles by the jet survived probability in and out of plane, Fig. 01 The parameter 
V2 reaches 11-12% in mid-central events (centrality 30-35%) and nicely follows the 
trends observed in the experiments [ 1111 |5J- Our result disproves the assertion that 
jet quenching models can not explain the measured vi at high px [ I12| . It also 
means that measured vi at high px is not sensitive to hydro/collective flow. 

5) Combining statements 1 and 4 we come to the conclusion that at intermedi- 
ate momentum range, 2-4 GeV/c, particles from early corona with the power law 
spectrum shape and with the smaller value of V2 dilute the exponentially falling 
spectrum from the core with hydro-like V2 value. This is why v-i at momenta below 
2 GeV/c rises with close to hydro calculation values, saturates at some value, and 
then starts to fall down [E3 to the values we present in Fig.|3| Would we develope 
the method to subtract corona contribution from the observed v-i value, the agree- 
ment with hydro calculation will be even better. At intermediate px region mesons 
are much more suppressed than baryons, this is why mesons azimuthal anisotropy 
goes down earlier than for baryons and, so called, constituent quarks t>2 scaling [ 
ITU] could be accidental. 

6) Undistorted production of high px particles from the corona region explains the 
lack of change of hadron distributions in forward jets with centrality and orientation 
relative to the reaction plane [ I17j . 

7) The same is true for di-jets. The yield of high px di-jets will be suppresed, lead- 
ing dominantly to the tangential di-jet prodiction with no change of jet properties. 
STAR collaboration published the results on Iaa in Au+Au collisions - the ratio 
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of away-side yield per trigger particle to the similar value from p+p collisions [I17j. 
We calculate Iaa within our model with T = 2.3/m/c for different centralitics. The 
increase of corona thickness by 2.3 fm keeps the value of Iaa at the level close to 
the experimental data, Fig. 0] 
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Fig. 4. Calculated Iaa for Au+Au, the curve, versus the number of participant 
nucleons, N par t- The width (sigma) of away-side jet is 0.5 radians. The cicles are 
experimental data from [I17j. 



8) In Fig. [3] we present PHENIX preliminary results on non-photonic electrons 
at high momenta shown at Quark Matter 2005 [I18[. We also plot the lower limit of 
Raa determined by particle production from the corona region during the formation 
time, solid line. Experimental points sit exactly on the line. These electrons origi- 
nated from D-mesons, charmed baryons and, probably, from bottom quark decays. 
The conclusion of Fig. [3] would be that at high p t the charm quarks are completely 
absorbed by dense medium, similar to the light quarks. 

9) In the case of a strong charm quark absorption, following the argument in our 
previous part 4, one can expect a significant V2 for non-photonic electrons. Experi- 
ment indeed found non-zero V2 [ HHj . 

10) Many interesting features may come from 3/ip production in nuceus-nucleus 
collisions. In our consideration we devide nucleus-nucleus collision in to two stages: 
before plasma formation and after. So, it should be different 3 ftps absorbtion pro- 
cesses at each stage of the collision. If in some momentum range 3 ftps are absorbed 
by the core, there should be non-zero v 2 too. 

11) Few speculative conclusions could be made from the beam energy dependence 
of T. At smaller beam energy the measured formation time rises, at 200 GeV T=2.3 
fm/c, at 62 GeV T=3.5 fm/c. At lower energy T should be even larger, such that 
sQGP can't be formed at all because of fast longitudinal expansion. Indeed, at 
20 GeV, see talk by C. Klein-Bosing at this workshop, Raa—^- Such behaviour 
goes in favor of the existence of strongly coupled 3- and 4-body bound states in 
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Fig. 5. Raa for non-photonic electrons in Au+Au versus the number of partici- 
pant nucleons, N part . The PHENIX preliminary experimental data are from [I18|: 
circles and squares are for p t =4.25 GeV/c and 4.75 GeV/c, respectively. Shown 
are statistical errors only. Solid line - the result of our calculation which puts the 
geometrical limit in case of strong charm quark absorption. 



sQGP phase, as it was proposed by Edward Shuryak and Ismail Zahed [J2J. Time, 
necessary to form such states, should be proportional to the mean distance between 
interaction points with color exchange. This distance, itself, should be inversally 
proportional to the square root of density of such interactions. If relativistic rise 
of the total nucleon-nucleon cross section comes purely from the parton hard scat- 
terings, we can estimate the relative contribution of hard scatterings in the total 
number of nucleon-nucleon collisions. At 20 GeV the NN total cross section is at its 
minimum - there is no hard scattering. At 62 GeV cross section rises by 5 millibarns, 
at 200 GeV - by 13 mb, at 5500 GeV - by 49 mb. Formation time of bound states 
should be proportional to one over the square root of these numbers. If we measure 
T=2.3 fm/c at 200 GeV, from rise of the NN total cross section, we estimate T=3.6 
fm/c at 62 GeV, T=1.2 fm/c at 5500 GeV, LHC energy. 

12) The existence of formation is a direct sign that sQGP is actually formed at 
RHIC. 

We present an orthogonal point of view on some experimental data. This view 
could be useful for better understanding of properties of a new matter produced at 
RHIC. Perhapes, it will bring clear vision or some new ideas on formation process 
and properties of a quark-gluon phase. 
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